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ABSTRACT 



The present paper reports on the RXTE observations of the binary X-ray 
pulsar 4U 0115+63, covering an outburst in 1999 March-April with 44 pointings. 
The 3—30 keV PCA spectra and the 15—50 keV HEXTE spectra were analyzed 
jointly for the cyclotron resonance features. When the 3—50 keV luminosity at 
an assumed distance of 7 kpc was in the range (5— 13)xl0 37 erg s -1 , harmonic 
double cyclotron features were observed in absorption at ~11 and ~22 keV, as 
was measured previously during typical outbursts. As the luminosity decreased 
below ~5xl0 37 erg s -1 , the second resonance disappeared, and the fundamental 
resonance energy gradually increased, up to ~16 keV at 0.16xl0 37 erg s -1 . These 
results reconfirm the report by Mihara et al. (2004) using Ginga, who observed a 
single absorption at ~16 keV in a minor (~ 10 37 erg s _1 ) outburst of this object. 
The luminosity-dependent cyclotron resonance energy may be understood as a 
result of a decrease in the accretion column height, in response to a decrease in 
the mass accretion rate. 

Subject headings: pulsars: binary, magnetic fields, cyclotron — X-ray: pulsars 
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1. Introduction 

Accreting binary pulsars are considered to have strong surface magnetic fields in the 
range of several times 10 12 G. One of the methods to accurately measure their fields is 
to observe cyclotron resonance scattering features (CRSFs) in their X-ray spectra, because 
the resonance energy E a and the magnetic field strength B are related with each other as 
E a = 11.6 -812(1 + Zg)^ 1 [keV], where B\ 2 is the magnetic field strength in units of 10 12 
Gauss and z g is the gravitational redshift. So far, CRSFs have been detected from ~ 15 X- 
ray pulsars, all in absorption, with balloons (e.g., Triimper et al. 1978), HEAO-1 (Wheaton 
et al. 1979; White et al. 1983), Ginga (e.g., Clark et al. 1994; Mihara 1995; Makishima et 
al. 1999), RXTE (e.g., Coburn et al. 2002; Heindl et al. 2001), BeppoSAX (e.g., Santangelo 
et al. 1999; Orlandini et al. 1998), and other missions. 

The recurrent transient 4U 0115+63, with the 3.6 s pulsations (Rose et al. 1979), is one 
of the X-ray pulsars whose CRSF has been studied in great detail. The optical companion is 
an 09e star, V635 Cassiopeae (Unger et al. 1998), with the orbital period of 24.3 days. The 
distance to 4U 0115+63 is estimated as 7 kpc (Negueruela and Okazaki 2001). Its CRSF was 
first discovered at ~23 keV in the HEAO-1 A4 spectra by Wheaton et al. (1979). Using the 
HEAO-1 A2 data obtained in the same outburst, White et al. (1983) suggested that the 23 
keV feature is in fact the second harmonic resonance, with the fundamental resonance at ~11 
keV. The suggested double harmonic structure was reconfirmed with Ginga by Nagase et al. 
(1991). Moreover, the third harmonic feature was found with RXTE by Heindl et al. (1999), 
and the fourth harmonic with BeppoSAX by Santangelo et al. (1999). In the BeppoSAX 
data, the resonance energies were 12.7, 24.2, 35.7 and 49.5 keV, or nearly harmonic ratios of 
1:1.9:2.8:3.9. Thus, 4U 0115+63 is one of the most suitable objects to study the physics of 
cyclotron resonance in the polar caps of binary X-ray pulsars. 

This source has been observed at various X-ray luminosity levels, and the fundamental 
resonance energy has been measured repeatedly at ~11 keV. However, in a minor outburst 
in April 1991 observed with Ginga, when the 2—60 keV luminosity (1.9 xlO 37 erg s _1 ) was 
~7 times lower than those typical at the outburst peaks, a drastic change in the CRSF was 
detected: instead of the familiar double absorption features at ~11 keV and ~22 keV, it 
exhibited a single deep and wide absorption at ~16 keV (Mihara 1995; Mihara et al. 1998; 
Makishima et al. 1999; Mihara et al. 2004, hereafter Paper 1). Such a large change in E a , 
presumably depending on the luminosity, had been observed neither from 4U 0115+63 itself 
previously, nor from other sources. 

In Paper 1, we tentatively concluded that the fundamental resonance energy, at ~11 
keV in normal outbursts, increased to ~16 keV, because a lower luminosity would make the 
accretion column shorter, and hence increase the magnetic field intensity at the column top. 
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However, there has remained an alternative possibility that the second harmonic resonance, 
normally at ~22 keV, decreased to ~16 keV in the 1991 outburst. Furthermore, even if 
the former interpretation is correct, the observed change in E a was considerably larger than 
is predicted by an accretion column model by Burnard et al. (1991). To make the results 
in Paper 1 less ambiguous, we need to more densely sample the spectra as the luminosity 
changes over a wide range. 

In this paper, we analyzed the RXTE data of 4U 0115+63 acquired in the 1999 March- 
April outburst, and studied the resonance energy as a continuous function of the luminosity in 
the rising and declining phases of the outburst. We have indeed detected the clear luminosity- 
dependent changes in the CRSF, and successfully reconfirmed the inference made in Paper 
1. All of the errors appeared in this paper are 90 % confidence levels. 

2. Observations 

The outburst of 4U 0115+63 to be utilized in the present paper was detected in 1999 
March with the All Sky Monitor (ASM) on board RXTE. Figure 1 shows the light curve 
of the whole outburst. The peak intensity, ~30 c s -1 , corresponds to approximately 400 
mCrab, which is typical of regular outbursts of this source. The BeppoSAX observation of 
this outburst led to the discovery of the four harmonic CRSFs by Santangelo et al. (1999). 

During this outburst, 44 short (~ 1 ksec) pointing observations were made with the 
Proportional Counter Array (PCA; Jahoda et al. 1996) and the High-Energy X-ray Timing 
Experiment (HEXTE; Rothschild et al. 1998) on board RXTE, starting on 1999 March 3 and 
ending on April 20. Heindl et al. (1999) utilized these datasets, and successfully detected up 
to three CRSFs. The PCA consists of five xenon proportional counter units (PCUs) with 
an energy range of 2—60 keV and a total effective area of ~7000 cm 2 . In order to extend 
the detector life time, some of the PCUs were not operated in the present observations. In 
this paper, we utilize all PCUs operated at each observation, as summarized in Table 1. 
The HEXTE consists of two arrays (cluster A and B) of four Nal/Csl scintillation counters, 
with an energy range of 15—250 keV and a total effective area of ~1600 cm 2 . However, one 
pulse height analyzer of cluster B failed in the early mission phase, so its effective area is 
approximately 3/4 of the nominal value. 

From the 44 pointing observations of the 1999 outburst in the RXTE archive, we selected 
34 data with low electron contamination; these are listed in Table 1. In this selection, we set 
the electron-rate threshold at 0.12 counts s" 1 (which is a typical value when the signal X-ray 
flux is high), in order not to discard good data by mistaking signal X-rays for electrons. In 
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order to complement the rather high threshold, we visually inspected the electron-rate light 
curves, and confirmed that none of the 34 datasets are contaminated by sporadic electron 
events. We then selected the Standard-2 mode data which have 16 s time resolution and 
129 energy channels, and used only the top layer of the PCA which has a low background. 
We employed FTOOLS v5.2 for analysis, and used the calibration data file epoch 3 and 4 
to make the energy response matrix. The background subtraction utilized so-called bright 
background model for all data. We have selected those data acquired when the offset angle 
of the source to the field of view center was < 0°.02, the earth elevation angle was > 10°, and 
the spacecraft was not within 30 minutes of an entrance from the South Atlantic Anomaly. 
The net exposure of each datasets, which was obtained after these screenings, is also given 
in Table 1. 

Since different observations used different sets of PCUs (Table 1), it is important to 
evaluate systematic errors associated with the PCA responses (particularly among different 
PCUs). Accordingly, we followed Wilms et al. (1999) and Coburn et al. (2001), and analyzed 
the Crab nebula data obtained on 1997 December 12 (for epoch 3) for 2.8 ~ 5.4 ksec, and 
1999 December 18 (for epoch 4) for 1.9 ~ 8 ksec, with the exposure depending on PCUs. By 
fitting these data with a single power-law of photon index 2.1, which is known to be a good 
approximation of the Crab spectrum in the PCA range, we examined how the data-to-model 
ratio at each energy bin scatters among the PCUs. Form this, we have set the systematic 
error of the PCU to be 1 %. The subsequent results do not change if we instead use 0.5 %. 



3. Analysis and Results 

3.1. Ratios to the Crab spectrum 

For each dataset listed in Table 1, we accumulated 3—30 keV PCU data into a single 
spectrum, and 15—50 keV HEXTE data into another. The background was subtracted in 
the way described in section 2. In order to grasp the spectral shape in a model-independent 
manner without being hampered by the instrumental response, we surveyed these spectra 
normalized to that of the Crab spectrum. The reference Crab data are the same as utilized 
to estimate the systematic error. The results are shown in Figure 2, where errors which 
propagated from the Crab spectrum are negligible. 

In Figure 2, the first data set acquired on March 3 gives a relatively featureless spectral 
ratio. When the flux increased and exceeded 710 c s _1 (PCU) -1 on March 5, the Crab ratio 
started to reveal two dip features at ~12 keV and ~25 keV. These can be interpreted as the 
fundamental and second harmonic, as have been observed repeatedly. The double absorption 
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feature persisted throughout the flare peak until March 29. However, when the flux decreased 
to <400 c s _1 (PCU) -1 on April 2, it turned into an apparently single absorption feature. 
Thus, the present data appears to reconfirm the Ginga results reported in Paper 1, that the 
CRSF appears in double absorption when the X-ray intensity is high, while it changes into a 
single absorption as the source becomes faint. In the RXTE data, the threshold is suggested 
to lie between 400 and 710 c s _1 (PCU) -1 . Below, we attempt to quantitatively confirm 
these inspections. 



It is known that typical continuum spectra of accreting binary pulsars can be approxi- 
mated with a power-law times exponential cutoff model. In this paper, we employ its updated 
version called NPEX (Negative and Positive power-laws with Exponential) model (Mihara 
1995; Makishima et al. 1999), and fit the background-subtracted spectra in comparison with 
the results reported in Paper 1. The NPEX model is written as 



where E is the X-ray energy in units of keV, A\ and a\ are the normalization and photon 
index of the negative power-law, respectively, Ai and 0:2 are those of the positive power-law, 
and kT represents the cutoff-energy in units of keV. In this paper, we tentatively fix a 2 at 
2.0, so that the positive power-law describes a Wien peak (Rybicki & Lightman 1979). 

Figure 3 shows the pulse-phase-averaged PCA+HEXTE spectra obtained on March 27 
and April 8, when the Crab ratio suggests the double and single CRSF, respectively (see 
Figure 2). We first attempted to fit the spectra with the NPEX model. In order to take 
into account possible over- or under-subtraction of background, we allowed the background 
normalization to vary (as described in the RXTE cook book), so as to minimize the fit chi- 
squared. We found the optimum normalization factor to be close to the nominal value of 1.0 
within a few percent. Even allowing this correction, the model left significant structures in 
data-to-model ratios as shown in Figure 3d and g. As a result, the fits remained unacceptable 
with xl ~ 62 and xt ~ 30 for the March 27 and April 8 spectra, respectively. Indeed, the 
March 27 spectrum exhibits two negative deviations at ~12 and ~22 keV from the NPEX fit, 
while that of April 8 shows only one negative feature around 16 keV. These results reconfirm 
the inference from the Crab ratios. 

In order to better reproduce the spectra, we next introduced a cyclotron absorption 



3.2. Analysis of representative spectra 
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(CYAB) factor which is given as 




r (WE/E a ) 2 



(2) 



(E - E a f + W 2 



where E a is the resonance energy, W is the width of the absorption structure, and r is 
the depth of the resonance (Mihara et al. 1990). We fitted the spectrum on March 27 with 
the NPEX continuum multiplied by two CYAB factors (hereafter NPEXxCYAB2 model), in 
which all the NPEX and CYAB parameters are left free except that the second CRSF energy 
is fixed at twice the fundamental energy. For the spectrum on April 8, we applied the NPEX 
multiplied by a single CYAB factor (hereafter NPEX x CYAB model). The PCA+HEXTE 
data were fitted simultaneously with the same parameters, but using another free parameter 
to adjust relative normalizations of the two instruments. The background normalization was 
again allowed to vary, by up to ±1.5 %. We did not include low energy absorption nor Fe-K 
line, since our spectra exhibit neither low-energy turn-offs nor an excess around 6.4 keV. 
The lack of significant Fe-K line is often noticed among transient X-ray pulsars with Be-type 
primary stars (e.g., Nagase 1989); in such a system, the matter to be accreted by the pulsar 
presumably forms a disk-shaped envelope around the Be star, and hence its solid angle as 
seen from the pulsar is much smaller than the case where the matter is captured from a more 
isotropic stellar wind. 

As shown in Figure 3e and f, the NPEX x CYAB model has successfully reproduced 
the PCA+HEXTE data on April 8, yielding reduced chi-squared of ~1.2. The background 
correction factor turned out to be 0.3 %. In contrast, the March 27 spectrum was not 
reproduced successfully (xt ~ 2.0) even by the NPEXxCYAB2 model, because of negative 
residuals around 35 keV (Figure 3c). This feature was removed by multiplying a third 
CYAB factor, of which the centroid energy is 32.8+1J keV. Since this energy is close to 
3 x E a i, the feature can be identified with the third harmonic resonance detected from the 
same outburst (Heindl et al. 1999; Santangelo et al. 1999). We hence fitted the March 
27 spectrum with the NPEX continuum multiplied by three CYAB factors, of which the 
resonance energies were constrained to have harmonic ratios (1:2:3). This model, hereafter 
called NPEXxCYAB3 model, successfully reproduced the March 27 spectrum with xl ~ 1-2. 
The derived fundamental resonance energy on March 27 is 10.4 ± 0.1 keV at an X-ray 
luminosity of 7.6 xlO 37 erg s" 1 , while that of April 8 is 14.2 ± 0.4 keV at 3.0 xlO 37 erg s^ 1 . 
The resonance energy has thus increased by a factor of 1.4 as the luminosity decreased by a 
factor of 2.5, although there still remains a possibility that the single CRSF at 14.2 keV on 
April 8 is in reality the second harmonic. 

While we tentatively fixed a 2 at 2.0, Santangelo et al. (1999) obtained a 2 = 0.41 ±0.05 
by analyzing the BeppoSAX data of the same outburst. This is probably due to difference 
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in the energy band used; we fitted the spectra in the 3 — 50 keV range, while Santangelo et 
al. (1999) in 9 — 100 keV. To examine this difference, and to examine whether the constraint 
of «2 — 2.0 is justifiable, we fitted the NPEX x CYAB3 model, with a 2 left free, to the 
PCA+HEXTE data on March 19a, which is one of those with the highest statistics. Then, 
this parameter has been constrained as a 2 = 1-60 ± 0.99, which includes the fiducial value 
of 2.0. More specifically, fixing a 2 at 2.0, at our best fit value of 1.60, and at the BeppoSAX 
value of 0.41, gave a reduced chi-squared of 1.29 {y = 78) , 1.33 {y = 77) , and 1.67 [y = 78), 
respectively. Based on these results, we hereafter retain our initial assumption of a 2 = 2.0. 

Apart from the continuum modeling, we must examine whether the CRSF parameters 
we have derived are consistent with the previous analyses of the same outburst (Heindl 
et al. 1999; Santangelo et al. 1999). We therefore analyzed the PCA+HEXTE spectra on 
March 11.36—11.40 (Obs. 7; Table 1), in comparison with those reported by Heindl et 
al. (1999) for a particular pulse phase (0.70—0.76 in terms of their phase convention) on 
March 11.87—12.32. We found some difference in the resonance energy; 10.61^ keV from 
our fit, while 12.401q35 keV by Heindl et al. (1999). The difference is likely to arise from 
the different modeling of the CRSF employed in these studies, because Heindl et al. (1999) 
used a Gaussian absorption model whereas we are using the CYAB factor. Accordingly, 
we re-fitted the same spectra with an NPEX model multiplied by two Gaussian absorption 
models. The fundamental cyclotron resonance was then obtained at E al = 12.3 ±0.1 keV, 
which is consistent with Heindl et al. (1999). The higher accuracy of our result is probably 
because we analyzed the pulse phase-averaged spectra. Other resonance parameters turned 
out to be consistent between the two works. We thus conclude that our results are consistent 
with those of Heindl et al. (1999) using the longer exposure, and have comparable accuracy. 



3.3. Analysis of date-sorted spectra 

For the moment, we mainly concentrate on the study of the 3—30 keV PCA spectra, 
because this energy range is most relevant to the CRSFs. In order to more continuously 
resolve the drastic spectral change between March 27 to April 8, we applied the NPEX model 
(without incorporating CYAB factors) to the daily-averaged PCU0+PCU2+PCU3 spectra 
over this period. Figure 4a shows ratios of the data to the best-fit NPEX model, which 
allow us to grasp the spectral changes objectively. The PCA background normalizations 
were corrected in the way as described in section 3.2. Since the ratio on March 27 and 29 
shows almost the same shape as Figure 3d, we may consider that the two CRSF persisted 
at least until March 29. On March 31, the second CRSF became less clear, and on April 2 
onward, the ratio reveals only a single absorption. In addition, the fundamental resonance 
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initially observed at ~ 11 keV moved to higher energies, up to ~ 15 keV, while there is no 
opposite trend such as the two features drifting toward lower energies. 

In order to quantify the results of these visual inspections, we have fitted the same 
set of spectra with the NPEXxCYAB2 (or NPEXxCYAB) model. The obtained best-fit 
parameters are summarized in Table 2, and the ratio of the data to the model are plotted in 
Figure 4b. The spectra from March 27 through April 2 all required the two CYAB factors, 
because the double CYAB fit gave a significantly better reduced chi-square (0.6~1.2) than 
the single CYAB case (>2). On April 4, the spectrum is roughly reproduced by the single 
CYAB model, although there remains a hump at ~23 keV which does not disappear by 
including the second CYAB factor. The successful single CYAB fit continued to the end, on 
April 8. These results quantitatively confirm our inference made above, that the fundamental 
CRSF energy increased from 10 to 15 keV over this ~10 day period as assumed in Paper 1, 
although details of the change are not yet resolved clearly. 



3.4. Analysis of intensity-sorted spectra 

Although the 3—30 keV spectrum apparently depends on the luminosity, the source 
exhibits significant intra-day variations, as shown in Figure 5. As a result, each spectrum 
in Figure 4, in reality, is an average over a relatively wide luminosity range. In order to see 
the luminosity dependence more clearly, we sorted the PCA data from March 27 to April 
8 into 8 intensity intervals, in reference to Figure 5. We again used data from PCU 0, 2 
and 3, which worked throughout this period. Then, we have repeated the same analysis as 
performed in section 3.3. 

Figure 6a shows the ratios of these intensity-sorted spectra to their best-fit NPEX 
models, to be compared with Figure 4a. The two CRSFs, at ~11 keV and ~22 keV, are thus 
observed clearly in the higher four intensity levels, fl through f4. As the intensity decreased, 
the second CRSF gradually became shallower. Finally, the second CRSF disappeared at 
level f6, and the fundamental CRSF started to move from ~11 to ~16 keV over levels f5 
through f8. From these results, we confirm that the ~16 keV single structure results from 
an upward shift of the ~10 keV fundamental CRSF. The spectral changes revealed here are 
consistent with, but clearer than, those seen between April 2 and April 4 in the date-sorted 
spectra. 

Using the NPEXxCYAB2 (or NPEXxCYAB) model, we quantified the CRSF param- 
eters as a function of the intensity. The fitting results are summarized in Table 3, and the 
ratios of the data to the model are shown in Figure 6b. The full fitting results are displayed 
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in Figure 7. The two CYAB factors have been required by the fl to f4 spectra. The sec- 
ond CRSF is not clearly visible in the f5 ratio (Figure 6a and Figure 7), but a large xt 
(~2.0) was obtained by the single CYAB fitting. We therefore applied the NPEXxCYAB2 
model to the f5 spectrum, and obtained a fully acceptable fit (Table 3). This means that 
the second harmonic is still present in the f5 spectrum, with a reduced depth (t2 — 0.2). 
From fl through f5, both r 2 and Wi decreased; that is, the second harmonic feature became 
narrower and shallower. Although this behavior is absent in the day-sorted results (Table 2), 
the difference can be attributed to the fact that each day-average spectrum is a mixture of 
different spectra corresponding to different intensities. The f6 through f8 spectra have been 
fitted successfully by the NPEXxCYAB model. 

These results unambiguously show that the center energy of the fundamental CRSF 
moved from 10 to 15 keV as the luminosity decreased. Thus, we can conclude that the 
16 keV spectral structure of 4U 0115+63 observed by Ginga (1991) and this work is the 
fundamental CRSF, rather than the second harmonic resonance which moved toward lower 
energies. The threshold between the single and double CRSF structures is found at a 3 — 30 
keV luminosity of ~4.2xl0 37 erg s _1 . 

3.5. Analysis of all PC A and HEXTE spectra 

Now that we have understood the basic behavior of the CRSFs, we now proceed to 
the analysis of the entire PCA and HEXTE data prepared in section 3.1. Since the Crab 
ratios on March 3 and 4 do not show clear CRSFs (Figure 2), we first attempted an NPEX 
fitting, but the model was unacceptable (xt — 20), because of the negative deviation at 
~ 15 keV. Then we applied the NPEXxCYAB model to these data, but the fit was still 
unacceptable (xt — I-?)- Finally, we fitted these data with the NPEXxCYAB2 model, and 
obtained reasonable fit with xt < 1-3- These results imply that the double CRSFs were 
already present in these early data sets, with the fundamental resonance at ~ 12 keV as 
given in Table 4. Thus, the data from the beginning through April 2 were generally well 
described with the NPEXxCYAB2 model. By jointly analyzing the PCA and HEXTE data, 
the second CRSF parameters have been significantly better constrained than using the PCA 
data alone. The derived parameters are listed in Table 4. 

In some fits, however, we were still left with rather large values of xt- F° r example, the 
NPEXxCYAB2 fits are rather poor (xt > 1.5) on March 15, 19a, 21a, 22, 27, 29 and 31 
(Table 4), often with negative residuals at ~35 keV. Like in Figure 3b, these residuals can 
be mostly removed by multiplying the model with a third CYAB factor. Typical parameters 
of the third CRSF, e.g., observed on March 19a, are E a3 = 35.2 ± 2.2 keV, W 3 = 9.3 ± 2.5 
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keV and r 3 = 0.6 ± 0.1. We can hence identify it with the third harmonic, after section 
3.2. The third harmonic persisted in the spectra for about half a month from March 15, 
with relatively constant parameters, till March 31 when it became unconstrained presumably 
due to insufficient data statistics in higher energies. Meantime, the second CRSF stayed at 
around 22 keV. 

The data from April 3 to the end did not show the higher harmonics, with the NPEXxCYAB 
model giving acceptable fits (< 1.4). In order to estimate the upper limits on the second 
harmonic feature, we attempted to fit the data obtained from April 3 to April 8 with the 
NPEXxCYAB2 model. In addition to the constraint of E a 2 = 2E a i employed so far, the 
second CRSF width W2 was tied to Wi, because W2 is close to Wi (except for April 1 and 2) 
when the spectrum exhibits the double features. The obtained upper limits on r 2 are given 
in Table 4. 

From these studies using the PCA and HEXTE data, we reconfirmed that the cyclotron 
resonance energies increase as the X-ray luminosity decreases. These results are consistent 
with those derived in section 3.4 from the intensity-sorted study of the 3 — 30 keV PCA 
spectra. Figure 8 summarizes the fundamental resonance energies, derived through both 
the date-sorted and intensity-sorted analyses, as a function of the calculated 3 — 50 keV 
luminosity at 7 kpc. 

4. Discussion 

We have analyzed the 34 PCA+HEXTE datasets, covering the whole 1999 March- 
April outburst of 4U 0115+63 in which the 3 — 50 keV source luminosity changed over L x = 
(0.17-14) x 10 37 erg s" 1 . When L x > 7x 10 37 erg s _1 , we observed the familiar double CRSFs 
with E & \ — 11 keV. As L x decreased across a rather narrow range of (5 ± 2) x 10 37 erg s _1 , 
the second harmonic resonance disappeared, and the fundamental resonance energy increased 
from E al ~ 11 keV to ~ 16 keV. These results reconfirm Paper 1, and unambiguously reveal 
that the resonance energy increases as the source gets less luminous. 

4.1. The resonance energy variations 

Except for some hysteresis effects between the rising and decay phases of the outburst, 
E a i changes roughly as a single- valued function of L x . In particular, the intensity-sorted 
3—30 keV PCA spectra and the date-sorted PCA+HEXTE data imply consistent results. 
Furthermore, the Ginga results reported in Paper 1 fall on the same E a \ vs. L x relation. 
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Therefore, the phenomenon is inferred to have a good reproducibility as a function of the 
luminosity. 

The CRSF is known to depend also on the pulse phase, presumably because different 
scattering regions characterized by slightly different magnetic field strengths are observed 
as the neutron star rotates. When studying luminosity-related changes in the resonance 
energy, such pulse-phase dependent effects must be considered. In the case of 4U 0115+63, 
however, the phase dependent variation of E a i in the 1999 outburst was relative small, ~ 10% 
(Santangelo et al. 1999). Furthermore, the phase- resolved analysis performed in Paper 1, on 
the two different luminosity levels, did not affect the main results from the phase-averaged 
spectroscopy. We therefore concentrate, in the present paper, on the phase-averaged analysis. 

As mentioned in Paper 1, luminosity-dependent changes of the accretion column height 
provide the most likely explanation to our results, because the column is expected to become 
taller as the source luminosity increases (e.g., Burnard et al. 1991). Assuming that the 
CRSF is formed at a height h r above the neutron star surface in the accretion column, and 
the magnetic field strength there follows the dipole law, we then expect 



with R NS the radius of the neutron star. The factor (1 + z g ) describes the gravitational 
redshift, but below we assume z g = for simplicity. Then, the relative resonance height, 
K/Rns, can be related to the resonance energy as 



where E is the resonance energy to be observed on the neutron star surface. 

Substituting the observed value of E & i into Equation (4), we have calculated h v /R^s 
as a function of the X-ray luminosity. The results are shown in Figure 9. Since there is no 
a priori knowing of E , we employed two different values of E ; 18 keV and 20 keV, the 
former close to the observed maximum value of E al . Except for some differences between 
the two assumptions on E , the height h T thus increases in a rough proportion to the X-ray 
luminosity, up to ~ 7 x 10 37 erg s _1 where h T appears to saturate at ~ 0.2i? NS . In addition, 
h r /Rws ma y approach a finite value of 0.03 ~ 0.07, instead of zero, as L x decreases below 
~ 2 x 10 37 erg s -1 . 

The results presented in Figure 9 may be compared to a theoretical prediction by 
Burnard et al. (1991), who calculated the height of the accretion column h top as 



E al oc (R NS + h r ) 3 (1 + z g ) 



i 



(3) 




(4) 




(5) 
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Here, Lf^ d is the Eddington luminosity along the magnetic field which is identical to the 
conventional Eddington Luminosity for a 1.4M neutron star, L^ dd = 2.0 x 10 38 ergs s _1 , and 
H± is the ratio of the Thomson cross section to the Rosseland-averaged electron scattering 
cross section for radiation flows across the magnetic fields. The dashed line in Figure 9 
shows this prediction, assuming H± = 1.23 (Paper 1) and -Rns — 10 km. Except for the 
observed saturation toward the highest luminosity and < 2 x 10 37 erg s _1 , the value of h r 
measured at each luminosity level thus corresponds to ~ 70% of the predicted h top . This 
is quite reasonable, because the resonance energy is expected to sample the magnetic field 
strength which is measured at, or slightly below, the top of the column (i.e., h T < h top ). 

So far, we have assumed that the observed X-ray intensity changes are caused solely 
by actual variations of the intrinsic source luminosity. However, the intensity may also 
be affected by other extrinsic factors, such as partial obscuration by materials around the 
accretion disk. Such effects are known, e.g., in Her X-l (Mihara et al. 1991). In the present 
outburst of 4U 0115+63, Heindl et al. (1999) detected occasional appearance of "mHz QPO", 
with a period of ~ 500 sec and intensity changes up to ~ 40 %, and suggested that the 
phenomenon is possibly due to source obscuration by some ionized materials. If so, we 
would need a caution in interpreting the results from our intensity-sorted analysis. 

In order to investigate effects of the mHz QPO on our results, we inspected PCA light 
curves acquired over the period used for our intensity-sorted analysis (March 27 through 
April 8). As shown in the inset to Figure 10, a relatively clear QPO was found on March 
31, although the period is about 1000 sec instead of the 500 sec reported by Heindl et al. 
(1999). Then, we sliced the light curves as shown there, and produced a pair of PCA spectra 
corresponding to peaks and valleys of the QPO. Figure 10 presents these spectra and their 
ratio, as well as residuals when fitted individually with the NPEX model. The resonance 
energy thus clearly increases at the QPO valley. By fitting the spectra with the NPEX times 
CYAB2 model (with W 2 again fixed at Wi), we have constrained the resonance centroid 
as 10.4lg;? keV at the QPO peak with L x = 8.2 x 10 37 erg s _1 , and 11.8±g;§ keV at the 
QPO bottom with L x = 5.2 x 10 37 erg s~ x . Thus, the energy shift is statistically significant. 
Furthermore, these two data points line up in Figure 8 closely on the general L x vs. E ai 
correlation from the day-sorted and intensity-sorted analysis. We hence conclude that our 
basic results remain unaffected, and that the QPO, at least on this occasion, is likely to 
reflect real changes in the intrinsic luminosity rather than some obscuration effects. This 
inference is reinforced by the relatively flat spectral ratio, because the ratio should show 
an increased low-energy absorption if the obscuration were due to neutral material, or Fe-K 
edge feature if it were highly ionized. 
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4.2. Behavior of the other parameters 

So far, many authors (e.g. Paper 1 and Coburn et al. 2002) studied relations among 
the continuum and cyclotron line parameters. Our work provides a valuable opportunity 
to investigate how these parameters in a single system change in correlated ways, when 
the resonance energy varies. Among various correlations, a particularly interesting one is 
that between r and the W/E a ratio; in fact, Coburn et al. (2002) found a positive correla- 
tion between them over a large sample of phase-averaged spectra of X-ray pulsars. Later, 
Kreykenbohm et al. (2004) noticed that the same correlation holds for the phase-resolved 
spectra of GX 301-2, and argued that the correlation may be explained if the accretion 
column has a tall cylindrical shape rather than flat coin-shaped geometry. As presented in 
Figure 11, the fundamental and second harmonic widths of 4U 0115+63 from the present 
observations both depend positively on their respective resonance depths, in agreement with 
the correlation found by Coburn et al. (2002). These results, together with the argument 
by Kreykenbohm et al. (2004), strengthen the tall cylindrical column geometry which we 
invoked in § 4.1. For reference, we did not find particular correlations between E a and kT, 
or between E a and W. 

Figure 12 shows the spectral parameters as a function of the X-ray luminosity. There, 
the positive r vs. W/E a correlation of Figure 11 has been decomposed into luminosity 
dependent changes in W/E a (panel a) and r (panels c, d). Up to ~ 4 x 10 37 erg s _1 , r x thus 
stays relatively constant, with a hint of mild increase presumably due to the increased column 
density in the emission region. Over the same luminosity range, the fractional width W\/E & i 
increases more clearly. This would not be due to changes in the Doppler broadening (e.g. 
Meszaros 1992), since the NPEX kT parameter shown in Figure 12b, which is thought to 
approximate the electron temperature in the emission region (Mihara et al. 2004; Makishima 
et al. 1999), stays rather constant. Instead, the increase in the W\/E a i ratio could be a result 
of the luminosity-correlated elongation in the accretion column, which would cause a larger 
range of magnetic field intensities to participate in the resonance formation. 

Beyond ~ 4 x 10 37 erg s -1 , both T\ and Wi/E a i start decreasing clearly, while the 
2nd harmonic resonance develops rapidly both in depth (t2, Figure 12d) and relative width 
(W2/Ea2 in Figure 12a). In short, the CRSF makes a transition from the single feature at 
low luminosities to the harmonic double feature at higher luminosities. Since this occurs 
approximately over the luminosity range where the rapid change in E al takes place, we 
consider that the single-to-double transition of the CRSF has the same origin as the resonance 
energy shift. 

The behavior of the fundamental and 2nd harmonic parameters at the single-to-double 
transition may reflect basic differences in their elementary processes. The fundamental 
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resonance has a very large cross section, but it acts as scattering rather than absorption, 
because an electron which is excited by absorbing a photon of energy ~ E al will immediately 
return to the ground state by emitting a photon of nearly the same energy. The second 
harmonic resonance, though with a much smaller cross section, will in contrast act as pure 
absorption, because the electron excitation/deexcitation in this case occurs via absorption 
of a photon of energy ~ 2E a i and cascade emission of two photons of energies ~ E & \ each. 
The emitted photons will fill up the fundamental resonance, making it shallower (so-called 
two-photon effects; Alexander & Meszaros 1991). In the present case, the accretion column 
may be effectively thick to the fundamental resonance (ti > 1) essentially at all the observed 
luminosities, but presumably optically thin (r 2 < 1) to the second resonance photons when 
the luminosity is low and hence the accretion column is short. As the source luminosity 
increases, the column becomes taller and opaque to the second resonance, leading to the 
emergence of the second resonance feature. At the same time, the two-photon effect would 
reduce Ti, just as seen in Figure 12c. A similar effect was observed in the INTEGRAL 
spectrum of GX 301 - 2 by Okada et al. (2004), who reported that n of its ~ 35 keV CRSF 
decreased toward higher luminosities. 



4.3. The case of X0331+53 

Although 4U 0115+63 thus provides the first clear example of luminosity-dependent 
changes of the cyclotron resonance energy, there has emerged another promising case. This 
is the transient X-ray pulsar X0331+53 (V0332+53), from which a very prominent CRSF 
was detected with Ginga at E a = 28.5 ±0.5 keV together with a hint of the second harmonic 
(Makishima et al. 1990). This was observed in an outburst when the 2 — 60 keV luminosity 
was ~ 2 x 10 37 erg s _1 at an assumed distance of 3 kpc. A re-analysis of the same data using 
the NPEX continuum have revised the value slightly to E a = 27.2 ± 0.3 keV (Makishima et 
al. 1999). 

The object entered a bright outburst from 2004 November, becoming considerably more 
luminous than was observed with Ginga (Swank et al. 2004; Remillard 2004). Observations 
with INTEGRAL (Kreykenbohm et al. 2005) and RXTE (Coburn et al. 2004) have clearly 
reconfirmed the fundamental CRSF, and further revealed the second and third resonances. 
The fundamental resonance energy obtained with INTEGRAL is 24.9 ±0.1 keV (Kreyken- 
bohm et al. 2005) at a 3 — 50 keV luminosity of ~ 8 x 10 37 erg s _1 . This means that a factor 
4 increase in the luminosity (from the Ginga outburst to the present one) is accompanied 
by a ~ 10% decrease in E & \. As mentioned by Mihara et al. (2004), the same effect was 
already visible between two pointings with Ginga. For comparison, the resonance energy 
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of 4U 0115+63 changes by a somewhat larger amount (~ 20%) across the same factor 4 
luminosity range (Figure 8). The archival RXTE data of X0331+53 are currently being 
analyzed, and the results will be reported elsewhere (Nakajima et al. , in preparation). 

Now, we have two examples of luminosity dependent changes in the resonance energy. 
However, there is an apparent counter example, namely Her X-l. In a large amount of RXTE 
data covering its "main on" phase, the 2 — 30 keV luminosity (at 5.8 kpc) of Her X-l varied 
by a factor of 2 over (1.8 - 3.2) x 10 37 erg s -1 , but the CRSF stayed at ~ 40 keV (Gruber et 
al. 2001). One possible interpretation of these results on Her X-l is that its luminosity swing 
was not large enough to reach the critical range (cf. Figure 8) where E a starts changing 
significantly. Alternatively, the luminosity-related changes in the CRSF may depend on the 
object, for some reasons which are yet to be detailed. 

In conclusion, our study using the RXTE data has confirmed the inference made in 
Paper 1, that the cyclotron resonance energy of 4U 0115+63 increases as the X-ray luminosity 
decreases. While this provide a new tool with which we can diagnose the accretion column of 
strongly magnetized neutron stars, it remains yet to be confirmed whether the phenomenon 
is common among this type of objects. 
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Table 1. The log of RXTE observations of 4U 0115+63 in the 1999 March-April 

outburst. 



PCA HEXTE 



Obs. 


Date 


Start /End Time a 


PCU. 


Exposure 


Ratc b 


Exposure 


Rate c 




(1999) 


(UT) 


No. 


[ks] 


[c s- 1 PCU" 1 ] 


[ks] 


[cs- 1 ] 


1 


Mar 3 


03:35/03:49 


all 


0.57 


556 ± 1.1 


0.18 


48.8 ± 0.9 


2 


Mar 4 


13:08/13:41 


0,1,2 


1.01 


466 ± 1.0 


0.15 


37.2 ± 0.8 


3 


Mar 5 


20:35/20:57 


all 


1.28 


713 ± 1.4 


0.44 


70.3 ± 0.6 


4 


Mar 6 


20:47/20:56 


all 


0.54 


806 ± 1.6 


0.18 


82.8 ± 1.0 


5 


Mar 7 


20:40/20:55 


all 


0.84 


866 ± 1.7 


0.29 


85.3 ± 0.8 


6 


Mar 9 


06:22/06:36 


0,2,3,4 


0.81 


918 ± 1.8 


0.20 


93.3 ± 0.9 


7 


Mar 11 


08:41/09:40 


0,2 


1.60 


1020 ± 2.0 


0.40 


97.7 ± 0.6 


8 


Mar 13 


18:58/19:41 


0,2,3,4 


0.67 


1011 ± 2.0 


0.21 


104 ± 1.0 


9 


Mar 14 


06:44/06:56 


0,2,3 


0.68 


1055 ± 2.1 


0.27 


105 ± 1.0 


10 


Mar 15 


07:53/08:17 


0,2 


0.79 


994 ± 2.1 


0.28 


101 ± 1.0 


11 


Mar 16 


09:33/09:46 


0,2 


0.25 


997 ± 2.4 


0.09 


99.9 ± 1.4 


12 


Mar 18 


07:49/08:05 


0,2,3,4 


0.37 


949 ± 2.0 


0.13 


98.5 ± 1.2 


13 


Mar 19a 


04:26/06:13 


0,2,3 


3.45 


924 ± 1.8 


1.07 


92.7 ± 0.4 


14 


Mar 19b 


12:51/13:12 


0,2 


0.33 


876 ±2.1 


0.10 


96.0 ± 1.5 


15 


Mar 20 


03:o4/Uo:o ( 


0,2 


0.88 


857 ± 1.8 


0.22 


or o _l_ n n 
80.2 ± 0.9 


16 


Mar 21a 


06:27/08:33 


0,2,3,4 


4.13 


875 ± 1.7 


1.46 


91.2 ± 0.3 


17 


Mar 21b 


12:48/13:27 


0,2,3,4 


0.93 


852 ± 1.7 


0.29 


88.9 ± 0.9 


18 


Mar 22 


04:22/05:28 


0,2 


3.11 


788 ± 1.6 


1.06 


82.0 ± 0.4 


19 


Mar 27 


07:53/11:01 


all 


4.02 


564 ± 1.2 


1.42 


60.2 ± 0.4 


20 


Mar 29 


00:12/03:01 


0,1,2,3 


4.26 


517 ± 1.1 


1.40 


49.2 ± 0.3 


21 


Mar 31 


04:11/06:24 


0,2,3 


4.11 


456 ± 1.0 


1.41 


41.5 ± 0.3 


22 


Apr 1 


03:27/04:19 


0,1,2 


0.27 


437 ± 1.2 


0.09 


37.8 ± 1.1 


23 


Apr 2 


04:09/06:14 


all 


3.63 


399 ± 0.8 


1.16 


32.7 ± 0.4 


24 


Apr 3 


02:30/02:45 


0,2,3 


0.50 


385 ± 0.9 


0.16 


30.8 ± 1.0 


25 


Apr 4 


04:08/06:11 


all 


3.47 


345 ± 0.7 


1.16 


26.9 ± 0.4 


26 


Apr 5 


11:09/11:27 


0,2,3 


0.89 


309 ± 0.7 


0.33 


23.0 ± 0.6 


27 


Apr 6 


02:25/04:28 


all 


3.80 


287 ± 0.6 


1.25 


20.8 ± 0.3 


28 


Apr 8 


01:13/03:27 


0,2,3,4 


4.66 


254 ± 0.6 


1.56 


17.3 ± 0.3 


29 


Apr 10 


03:59/06:05 


all 


3.48 


210 ± 0.5 


1.18 


12.1 ± 0.3 


30 


Apr 12 


02:16/04:28 


0,1,2,3 


4.12 


146 ± 0.3 


1.36 


8.6 ± 0.3 


31 


Apr 14 


01:14/03:22 


all 


3.94 


111 ± 0.2 


1.31 


6.8 ± 0.3 


32 


Apr 16 


21:58/22:12 


0,1,2,4 


0.83 


60.1 ± 0.2 


0.28 


2.9 ± 0.6 


33 


Apr 18 


21:56/22:10 


0,2,3,4 


0.82 


24.9 ± 0.1 


0.28 


1.9 ± 0.6 


34 


Apr 20 


01:14/01:34 


all 


0.64 


12.4 ± 0.1 


0.22 


0.7 ± 0.6 



a Start and end time (UT) of the PCA observations. 
b In the 3-30 keV energy range. 

c Count rates of HEXTE cluster A in the 15-50 keV energy range. 



Table 2. Summary of the best-fit parameters of the NPEXxCYAB(xCYAB) model, determined by the date-sorted 

PCA spectra from March 27 thorough April 8. 



Date 


(ph s 


- 1 koV- 1 ) 




ai 




(ph s 


- 1 kcV" 1 ) 


kT 
(keV) 


(keV) 




Wi 
(keV) 




ti 




w 2 

(kcV) 


T"2 


L = 


xl 


Mar 27 





50 +0.45 
° -0.22 


1 


49+° 


53 
49 


7 


9 „+1.15 
a -1.83 


4 91+ 1 12 
-0.43 


10.6+° 


1 

1 


+ 0.17 
°' al -0.15 





89±° 


05 
05 


6.80±l;f 9 


1 12+° 51 


6.79 


0.62 


Mar 29 





- 7 +0.50 
°' -0.25 


1 


<*±°o 


53 
49 


7 


00 + O.82 
S8 -1.24 


4 61+ 061 
* ,D1 -0.28 


10.6+° 


1 
1 


6 59+ 27 
— 0.25 





**±°0 


06 
04 


5.9l±l;5 7 


91 + ° 34 
u - al -0.19 


6.04 


0.89 


Mar 31 





„ + 0.34 
oo -0.20 


1 


68l° 


40 
38 


6 


, + 0.21 
1 -0.22 


5.22(flxed) 


10.7+° 


2 
2 


_ 7 „ + 0.29 
7 ' 7S -0.26 





95+° 


03 
04 


7 4 7 +0.69 
0.67 


1 03+ 07 
luJ -0.08 


5.25 


0.75 


Apr 2 





„.+0.25 
D -0.22 


1 


wis 


28 
37 


5 


9K +0.13 

^"-o.is 


5.39(fixed) 


11 5+° 


3 
3 


10 25+ - 76 


1 


«tg 


05 
06 


6.33±i;gl 


59+ ' 15 
u -° -0.20 


4.61 


1.13 


Apr 3 





O1 -0.27 


1 


79+° 


53 
47 


6 


r,n + 0.90 
oa -1.10 


, OO + 0.70 


13.6+5 





8 94+ 2 ' 27 
° ,y -1.63 





88±° 


22 
17 






4.24 


0.90 


Apr 4 





^+0.37 
°'-0.21 


1 




41 

38 


5 


„„ + 0.45 
8U -0.49 


3 Q5+ 32 
j.»o_ 22 


13.8±g 


4 


9 05+ 1 ' 16 
M ' UO -0.98 





93 ;° 


11 
10 






3.85 


0.93 


Apr 5 





,,+0.22 


1 


«±° 


42 
38 


3 


no + 0.88 
a,:> -1.01 


4 96+ 2 ' 59 
-0.83 


13.2+; 


8 
6 


, , 9 „ + 2.89 
11.28_ 2 02 


1 


27+° 
z ' -0 


26 
24 






3.42 


0.60 


Apr 6 





., + 0.28 
M -0.16 


1 


73+° 
'°-0 


40 
37 


4 


89j:°;*° 


3.92t«'j 2 


13.71° 


4 


8 .98±i;S 





91±° 


11 
10 






3.19 


0.97 


Apr 8 







1 


14+° 


30 




-0.48 


Jy 7 — 0.29 


15.0+° 


3 
4 




1 


oe ~ 


16 
15 






2.78 


0.94 



a xlfJ -3 

The second resonance energy Ea2 was fixed at 2 X E a l- 
°10 37 erg s" 1 in 3-30 keV 

I— > 

I 



Table 3. The same as Table 2, but for the intensity-sorted PCA spectra. 



Region 


Ai 


c*l 




feT 




W X 


Tl 


w 2 


T2 








(ph s- 1 keV- 1 ) 




(phs-^koV- 1 ) 


(kcV) 


(kcV) 


(kcV) 




(keV) 









fl 





42+°0 


30 
.17 


1 


27+ 


45 
.43 


7 


.20±° 


39 
.46 


5. 76(fixcdJ 


10. 


4+0 
-0 


1 

2 


s 7Q+016 
°- ,a -0.14 





.9510 


03 
.03 


7 QQ+0.52 
'• au -0.49 


1 50+ 04 
1 ' ou -0.05 


7 


.47 


0. 


84 


f2 


0. 


58l»- 


45 
24 


1 


■6ltS 


47 
.46 


8 


.7610. 


28 
.34 


4.59(flxod) 


10. 


6+0 
D -0 


1 
1 


5 97+0 16 
°- a '-0.14 





.8510 


02 
02 


5 94+ 40 
° ,y -0.39 


o-97l° :° I 


6 


.81 


0. 


72 


f3 


0. 


62+0, 


49 
.27 


1 




48 
.47 


8 


.4610. 


25 
.30 


4.54(flxcd) 


10 


6+0 
D -0 


1 
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Table 4. The best-fit parameters of the date-sorted spectra, incorporating the HEXTE data. 
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Fig. 1.— (top) The 2-10 keV RXTE ASM light curve of 4U 0115+63 in the 1999 March- 
April outburst, (middle and bottom) The 3-30 keV PCA and the 15-50 keV HEXTE light 
curves obtained from the 33 pointing observations listed on Table 1. 
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Fig. 2.— Background-subtracted PCA and HEXTE spectra of 4U 0115+63 in the 1999 
outburst, shown normalized to those of the Crab Nebula measured by respective instruments. 
Above plot exclude the same observation date(March 19b, 21b), low exposure data (March 
16, 18 and April 1). 
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Fig. 6. — (a) The same as Figure 4a, but for the intensity-sorted spectra defined in Figure 
5. (b) The intensity-sorted PCA spectra, each normalized to the best-fit NPEXxCYAB2 
(fl-f5) or NPEXxCYAB (f6-f8) model. Data are presented in the same manner as Figure 4. 
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the CYAB factor (s). 
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Fig. 8. — The obtained fundamental cyclotron energies shown against the 3 — 50 keV lu- 
minosity. The open symbols represent the fundamental energy of the double absorption 
features, and the filled symbols the energy of the single CRSF. All the date-sorted data in 
the brightening (circles) and declining (triangles) are presented together with the intensity- 
sorted (squares) results. The Ginga results (Paper 1) are also included with asterisks. 
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Fig. 9. — The same as Figure 8, but the resonance energy is converted through equation 4 
into the height (h T ) at which the resonance occurs. The dashed line represents the value of 
the column top height, h top , calculated by equation 5. 
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Fig. 10. — Spectral changes associated with a 1 mHz QPO observed on March 31. The top 
panel shows the background-subtracted PCA spectra accumulated at peaks and bottoms of 
the QPO, as specified on the 2 — 30 keV light curve (inset). The middle panel is the residuals 
from NPEX fits to the two spectra. The bottom pannel show the QPO bottom-to-peak ratio. 
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Fig. 11. — The CRSF depth r against the fractional width W/E a . The plotted points are 
obtained from the date-sorted and intensity-sorted analyses. 
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Fig. 12. — The dependence of the CRSF and NPEX parameters on the 3 — 50 keV luminosity, 
obtained from the 34 date-sorted data and the 8 intensity-sorted data. The circles and tri- 
angles represent those spectra which are described by the double and single CRSF factor (s), 
respectively, (a) The W j Ea ratio for the fundamental and second harmonic resonances, (b) 
The NPEX kT. (c) The fundamental resonance depth, t±. (d) The second harmonic depth, 

T~2- 



